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Abstract 

Using an analytical parameterization for the behavior of the x slope of the 
structure function F2 at small x in perturbative QCD, at the leading twist 
approximation of the Wilson operator product expansion, and applying a flat 
initial condition in the DGLAP evolution equations, we found very good agree- 
ment with new precise deep inelastic scattering experimental data from HERA. 
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1 Introduction 



The measurement of the deep-inelastic scattering (DIS) structure function F2 |T|, 0, 
and the derivatives ciFs/d ln(g2) g |^ | and d\nF2/ d\n{l/ x) |, | in HERA have per- 
mitted the access to a very interesting kinematical range for testing the theoretical ideas 
on the behavior of quarks and gluons [^] carrying a very low fraction of momentum of the 
proton, the so-called small x region. In this limit one expects that nonperturbative effects 
may give essential contributions. However, the reasonable agreement between HERA data 
and the next-to-leading (NLO) approximation of perturbative QCD has been observed 
for > 2 GeV^ (see review Ref. and references therein) and, thus, perturbative 
QCD could describe the evolution of F2 and its derivatives up to very low values, 
traditionally explained by soft processes. It is of fundamental importance to find out the 
kinematical region where the well-established perturbative QCD formalism can be safely 
applied at small x. 

The standard program to study the x behavior of quarks and gluons is carried out 
by comparison of data with the numerical solution of the Dokshitzer-Gribov-Lipatov- 
Altarelli-Parisi (DGLAP)Q equations |^ by fitting the parameters of the x profile of par- 



tons at some initial Qq and the QCD energy scale A |T^, |Tl|. However, for analyzing 
exclusively the small x region, there is the alternative of doing a simpler analysis by using 
some of the existing analytical solutions of DGLAP in the small x limit |jl2|-[l^. This 



was done so in Ref. [0 where it was pointed out that the HERA small x data can 
be interpreted in terms of the so-called doubled asymptotic scaling (DAS) phenomenon 
related to the asymptotic behavior of the DGLAP evolution discovered many years ago 



0. 



The study of Ref. |T2[ was extended in Ref. [0, |T^, |T5| to include the finite parts 
of anomalous dimensions of Wilson operators and Wilson coefficients]^, which predicts 
1^ , |15|] the small x asymptotic form of parton distributions (PD) in the framework of the 



DGLAP equation starting at some Qo with the flat function: 

fa{Ql) = (hereafter a = g, 51), (1) 

where fa are the parton distributions multiplied by x and Aa are unknown parameters to 
be determined from data. 

/^From now on, we refer to the approach of Ref. [|T^ [T^, IT3 as generalized DAS 



approximation. In generalized DAS the flat initial conditions writen in Eq. (|l|) determine, 
as in the standard case |T^, the basic role of the anomalous dimensions singular parts 
while the contribution from finite parts of anomalous dimensions and also from Wilson 
coefficients can be considered as corrections which however are important to have better 



agreement with experimental data [0]. In the present work, similarly to Refs. |T2[-|T5 
the contribution from the non-singlet quark component was neglected. 

The usage of the flat initial condition given in Eq. (|I]) is supported by the actual 
experimental situation: low-Q^ data [|l^, |l], |18|, H are well described for < 0.4 GeV^ 

^ At small X there is a different approach based on the Balitsky-Fadin-Kuraev-Lipatov (BFKL) equa- 
tion whose application is out of the scope of this work. However, sometimes we use below the 
BFKL-based predictions in the discussions and for the comparison with our results in generalized DAS. 

^ In the standard DAS approximation |l6| only the singular parts of the anomalous dimensions are 
used 



1 



by Regge theory with Pomeron intercept ap{0) = \p + 1 = 1.08, closed to the standard 
(ap(0) = 1) one. The small rise observed in HERA data H [18|, |1^, |[ at low Q"^ can 
be naturally explained by including higher twist terms (see [jl5|, |2^). Moreover, HERA 
data 1^, |T8|, ^ with > 1 GeV^ are in good agreement with the predictions from 



the GRV parton densities |]TT| which supports our aim to develop an analytical form for 
the parton densities at small x because, at least conceptually, our method is very closed 
to the GRV approach. 



The purpose of this work is to extend the study of Ref. [0] to compare the predictions 
from the generalized DAS approach with the new precise HI data for the F2 x slope. 

The article is organized as follows. In Section 2 we address the present situation with 
experimental data for the slope d\nF2/d\n{l/x) and we shortly review some approaches 
to describe them. Sections 3 and 4 contain, for completness, a compilation of the basic 
formulae in the generalized DAS aproximation from Ref. [1^ needed for the present 



study. In Section 5 we compare our predictions for the derivative d\nF2/d\n{l/x) with 
experimental data and discuss the obtained results. 



2 The slope dliiF2/dlii{l/x): experimental data and 
QCD phenomenology 

Various groups have been able to fit the available data (mostly separating the low 
and high regions) using a hard input at small x, x~^^ A > 0. This is clearly different 
from the flat input in the DAS approach of Refs. [T^-[|l^, also describing reasonably 
well the experimental results. In some sense it is not very surprising because modern 
HERA data (at large Q"^) cannot distinguish between the behavior from a steep input 
parton parameterization at quite large Q^, and the steep form adquired after dynamical 
evolution from a flat initial condition at quite low values. 

Moreover, for the evolution based on the full set of anomalous dimensions obtained 
at X — in Ref. |2T| within the BFKL formalism 0, the results weakly depend on the 



form of the initial condition p2l , preserving the steep ones and changing the flat ones. In 
the case of working with anomalous dimensions at fixed order in a^, the initial conditions 
are important when the data are considered in a wide range and it is necessary to 
choose adecuately the form of the PD asymptotics at some Q^. 

As we have already discussed in the introduction, the usage of a flat initial condition 
leads to the (generalized) DAS approximation [|T^-[]T3| (see the compilation of the most 
important formulas in the following two sections). An alternative to this is the choice of 
a steep initial condition at some quite large Q^: fa{x,Ql) ~ x'^" (subindex c stays for 
constant), that leads to the following Q^-dependence for fa{x,Q'^) [^-|^ (if x'^" » 
Const): 

fa{x,Q') M+(1 + A„g2) 



faix,Ql) M+(1 + A„Q2)' 
where Ml^{l + Ac, Q^) is the analytical continuation (from integer n to real 1 + Ac values) 
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of the '+' component of the Melhn moment of fa{x,Q'^): 

Ma{n,Q^) = f'dxx'^-^faix^Q^) (3) 
Jo 



When >> Const, the slope Ac should be Q^-independent |2^, and the whole 
Q^-dependence for fa{x, Q^) comes from the factor ~ M^{1 + Ac, Q^) in front of x'^" in 
Eq. dl) . Approximations similar to Eq. have been successfully applied in studying 
the Q^-dependence of HERA data at large values (see Ref. ||2^ and references therein). 



Considering separately the low Q"^ region, it is also possible to have a good agreement 
between F2 data and its Regge-like behavior Q]. Indeed, F2 at Q"^ ^ can be extracted 
from the relation: 

^ -crir, (4) 



where is electromagnetic coupling constant and a^o/ is the total (virtual) photopro- 
duction cross-section. 

A large amount of experimental data on hadronic cross-sections for many different 
processes shows an universal rise at large energies that gives the possibility to parameterize 
all these cross-sections as the sum of two different components 



72 = Ap.--w-i + A^."«w-\ (5) 



0" 

being s the center of mass energy squared. The constants Ap and A^i are process- 
dependent magnitudes and the intercepts Q;p(0) ^ 1.08 and a^lO) ~ 0.5 (see [^) are 
universal process-independent constants. The first and second terms in Eq. (^ correspond 
to (soft) Pomeron and Reggeon exchange, respectively. 

/^From Eqs. (^) and (|^) one immediately obtains that for ^ 

F2(x,g^) ~ x"^ and, hence, fa{x,Q^) ~ x'^ (e = ap(0) - 1 ^ 0.08), 

because s = Q'^/x at small x. 

There were many attempts to study the whole region in the framework of Regge- 
asymptotics (see, for example, the reviews in Ref. 0). As the reports in Ref. contain 
quite numerous sets of models, we will restrict ourselves to discuss below only two of 
them. 



In Ref. I^Ol the fit to F2 experimental data has been done under the approach 

/,(x,Q2)^^-A(Q^) (6) 

and a fast changing of the A((5^) has been found in the transition range ~ 5-^10 
GeV^. Unfortunately, it is rather difficult to reconcile for the whole range the Regge- 
like behavior given by Eq. (^) with DGLAP evolution. Some progress along this line has 
been done in Ref. that was also based on the flat initial conditions as given by Eq. 
(0). However, in Ref. |27| the PD structure is limited by the Regge-like form of Eq. (^, 



permiting to recoincile it with DGLAP evolution only separately at low Q^, where X{Q^^ 
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is close to (or to e), and at large Q^, where X{Q^) ~ Ac- The structure function F2 and 
parton distributions have been obtained in Ref. ||2^ for the whole range only as a 
combination of these two representations. 

For other types of models (see |]32|, the phenomenological Q^-dependence of A(Q^) 



has been introduced in the form: 



\{Q') = e\l + 



with a fitted constant c, that produces the soft values of the slope A(Q^) close to e at low 
Q2 and the hard ones A(Q2) ~ ~ 0.2 ^ 0.3 at > 20 GeV^. 

Very recently new precise experimental data on A(Q^) has become available 0. The HI 
data points are shown in Fig. 1 where one can observe that for fixed Q^, A is independent 
on X within the experimental uncertainties, in the range x < 0.01. Indeed, HI data are 
well described by the power behavior [|5|: 

F,{x,Q') = Cx-'^^"\ (7) 

where A(g2) = d\n{Q^/A^) with C ^ 0.18, a ^ 0.048 and A = 292 MeV. The linear rise 
of A with InQ^ given by Eq. (0) is plotted in Fig 2. 

A similar analysis for extracting X{Q^) as a function of x has been carried out by the 
ZEUS Collaboration. As it is possible to see in Fig. 8 of Ref. 0], the ZEUS data for 
A((5^) are compatible with a constant value ~ 0.1 at < 0.6 GeV^, as it is expected 
under the assumption of single soft Pomeron exchange within the framework of Regge 
phenomenology. For the case of HI, this behavior can also be inferred from the new 



preliminary HI data |34] at quite low values of Q . 



Let us point out that, even though our results in the framework of the generalized 
DAS approximation (Eqs. (^)-(|TTp below) do not have an explicit power-like behavior, 
they really mimic a power law shape over a limited region of x and (see Section 4). 
In adition we have observed earlier [|l^] that the x dependence of the effective slopes in 
generalized DAS is not strong and the F2 effective slope was in good agreement with old 



(less precise) HI data |jl|. We repeat in Section 5 the analysis performed in Ref. ||Tj] but 
now using the new precise HI data for the slope 0. 

3 dependence of F2 and parton distributions in 

generalized DAS approximation 

To begin with, we shortly write down the results of the generalized DAS firstly pre- 
sented in Ref. [|T1|. The small x behavior of the parton densities and F2 at NLO has the 



form^ : 

faiz,Q') = f^iz,Q') + f -{z,Q') and 

^ Hereafter z = x/xq, where xq is a free parameter which hmits the apphcabihty range of formulae 
(PI) -(pi]) and can be fitted from experimental data together with the magnitudes of gluon and sea quark 
distributions at Qg. As it has been shown in Ref. the fits to F2 HERA data depend very slightly on 
the concrete xq value. 
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/-(z,g2) ~ expM_(l)s-Z}_(l)p) + 0{z) 
f+{z,Q') ~ /o(a)exp(-rf+(l)s-D+(l)p) + 0(p) 



(1 - 4_(l)a,(Q2))^l(^ + 20asiQ'] 



(9) 

+ 0{p) (10) 



where e = X]f e^// is the average charge square of / effective quarks, as = as/(47r) and 



In 



p = as{Ql) - as{Q^ 



D± = d±± - —d± 



a = 2y {d^s + D^p)\nz 



{d^s + D+p) 



a 



In z 



2\n{l/z) ' 



:i2) 



and /Sq and Pi are first two terms of QCD /^-function. 

The components of the leading order (LO) AD d„{n) and singular parts d^ and regular 
ones d^{n) of the LO AD d^{n) = d^/{n — l) + d^{n) have the following form (at n ^ 1): 



12 

To 



d+il) = 1 + 



20/ 



djl] 



16/ 



27/?o ' ' 27/?o 

The corresponding NLO components can be represented as follows: 

412 



(13) 



d^M] 



dlJl) 



dljl) 



f 



/3o 



36C3 + 33C2 



-20 
1643 



di_ 







12 9 



/ 



68 



13 



4(2 - — / 
9 ^ 243 ■ 



134 13, -,g 80, 

— -12C2-^/ , dlM) = ^f , 

4f2C3-3C2 + ^ + /f4C.-^ + i^/ 



9/3o 





d'^Jl] 



-3 1 



/ 



(14) 



Some interesting features of the results given in Eqs. (|8|)-(P^ are summarized below: 

• Both, the gluon and quark singlet densities given above are presented in terms of 
two components ('+' and '— ') which are obtained from the analytical Q^-dependent 
expressions of the corresponding ('+' and '— ') components of PD moments. 

• The '— ' component is constant at small x, whereas the '+' component grows at 

> Qo ~ exp(cr), where a contains the positive LO term |(i_|_|sln(l/2;) and the 
negative NLO one \D+\phi{l/z) (see Eq. (|12D). So, the most important part from 
the NLO corrections (i.e. the singular part at x — > 0) is taken in a proper way: 
it comes directly into the argument of the Bessel functions and does not spoil the 
applicability of perturbation theory at low x values. 
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4 dependence of the slope dXwFij dXwiXj x) in gen- 
eralized DAS approximation 

The behavior of the parton densities and the structure funcion F2 within generahzed 
DAS approach, given by Eqs. (|D-(|ll|), can be represented by a power law shape over a 
hmited region of x and [0, : 

/,(a;,Q2) ^ a;-Af^(-,Q^) and F^ix.Q'') ~ x""- (^'^'^ 

Because djdXwx = d/dlnz, the effective slopes can be obtained directly from Eqs. 
(§)-(|ll|). They have the form: 

m.Q^) - #^-p ^^^^^ 



f,{z,Q') + {2f)/3as{Q')-f,{z,Q') 



It is interesting to emphasize that from Eq. ([T5|) one obtains that the gluon effective 
slope Xg^f is larger than the quark slope A^-^-^ which is in excellent agreement with 
MRS H and GRV JTll analyses (see also Ref. Jig). 

On the other hand the effective slopes X^J^ and X^p2 in Eq. (plSf ) depend on the 
magnitudes of the initial PD, Aa, and also on the chosen input values Qg ^'^d A. However, 
at quite large Q^, where the "— " component is negligible, the dependence on the initial PD 
disappears, having in this case, for the asymptotic behavior, the following expressions^: 



h{a) ^ _41n(l/^) 
.eff,asr = h{o){l-d%_{l)as{Q')) + 2Qas{Q')h{a)/p 

' ^ h{a){l-dl_{l)as{Q^)) + 2Qas{Q^)h{a)/p 



Xl^f'^^{z,Q') + 



\n{l/z) 



where the symbol ~ denotes that an approximation was done in the expansion of the 
modified Bessel functions /^(cr) (n = 0, 1,2). These approximations are accurate only at 
large a values (i.e. at large and/or small x). 

^ The asymptotic formulae given in Eq. work quite well at any > Qq values, because at 

= Qq the values of A^-''-'^ and A^^j'^ are equal zero. The use of approximations in Eq. (|l^) instead of 



the exact results given in Eq. (15) underestimates (overestimates) only slightly the gluon (quark) slope 
at > Qo- For the F2 case, the similarity of X^l^ and A^^j'^'"'* values is shown in Fig 1. 



6 



Finally, we note that at LO the F2 slope A^^"'^'"'* is equal to the quark slope A^-^-^ 
it coincides, for very large values of a, with the result from Ref. [p6| 
input (see also the first article in Ref. 0). At NLO, X'^f'"-^ 
gluon slopes but closer to the former (see Fig. 3 in Ref. 



and 

in the case of a fiat 
^e/^/,ai, jjgg between the quark and 



Tl 



5 Comparison with experimental data 

Using the results of previous section we have analyzed HERA data for the slope 
(ilnF2/(iln(l/a;) at small x from the HI Collaboration^ 

Initially our results for \^p2 depends on five parameters, i.e. Qq, a;o, Ag, Ag and 
Aj^(f=4). In our previous article we have fixed Aj^|g-(f=4)=250 MeV which was a 
reasonable value extracted from the traditional (higher x) experiments. All other param- 
eters have been fitted and good agreement with F2 HERA data has been found for Qq ~ 1 
GeV^ (all results depend very slightly on xq). 

In this work we take Aj^|g-(f=4)=292 MeV in agreement with the more recent HI 
results IP and other analyses (see Ref. |^ and references therein) and we fit directly the 
slope (ilnF2/(iln(l/a:;) data using Eg . (pTSf) . The result is shown in Fig. 1. For comparison 
we also plot the curves from a fit to F2 data in Ref. where the value 250 MeV was 
used. The results are very similar and demonstrate the very important feature of an 
approximate x-independence of as given by Eq. ([T5|) , which is in agreement with HI 
data 0. 

We also present in Fig. 1 the asymptotic values for the slope Xpl'""^ as obtained from 
Eq. (p!6|). The agreement with data and with the other curves is also rather good if 
one takes in consideration that, in this case, there is no fit involved because the only 
parameters entering Eq. (|TBp are the fixed values Aj^(f=4)=292 MeV and xq = 1. 

Thus, the tiny x dependence of the slope X'p2{x, Q"^) in the considered region of x and 
support the possibility to successfully use our generalized DAS approximation in the 
type of X independent analysis for the F2 slope. 

Fig. 2 shows the experimental data for ^^^2 the corresponding HI parameteri- 
zation written above in Eq. (|^). We have also plotted the result from Eq. (p!6D and 



one from Eq. (0) using the parameters from our previous article [|14| as it was presented 
before in Fig. 1. For both cases, we give it for two representative x values. 

/^From the visual inspection of Fig. 1, one can notice that the boundaries and mean 
values of the experimental x ranges ^ increase proportionally with Q^, what is related to 
the kinematical restrictions in the HERA experiments: x ~ 10~'^-Q^ (see Refs. 0, 0, ^ 
and, for example. Fig. 1 of [|]). 

Fig. 3 shows HI experimental data for X'p2 and the HI parameterization (Eq. 
(0)), as given in Fig. 2, but this time in comparison with the asymptotic values Xp2''^^ 
calculated from Eq. (0) using x = a ■ 10""^ ■ with a = 0.1, 1 and 10. One has a 
reasonable agreement with HI data for > 2 GeV^ using a between 0.1 and 1 (the two 

^ In this work we only use HI data The prehminary ZEUS data for the slope (ilnF2/(iln(l/a::) 
are only available through points in Figs. 8 and 9 of Ref. They shown quite similar properties in 
comparison with HI data |^]. Unfortunately, the ZEUS numerical values are yet unavailable and we can 
not analyze them in the present work. 
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lower dashed curves in Fig. 3), which approximately corresponds to the middle points of 
the measured x range. 



6 Conclusions 



We have studied the dependence of the slope = d\nF2/d\rL{l/x) at small x in 
the framework of perturbative QCD. Our results are in good agreement with new precise 
experimental HI data |^ at > 2 GeV^, where perturbative theory can be applicable. 

Despite the fact that our approach, which can be called generalized doubled asymptotic 
scaling approximation, is based on pure perturbative grounds: a fiat initial conditions at 
Qq ~ 1 GeV^ and dynamical evolution to > Ql (it is very close conceptually to GRV 
approach but contains exact analytical evolution), it can be reasonable applied for the 
new precise data of the slope X^J-/ ■ 

The agreement between X'p2 data and perturbative QCD has been already observed 
by H10 and ZEUS [^] collaborations. The obtained linear rise of A(Q^) with InQ^ (see for 
example Figs. 2 and 3), parametrized by HI as in Eq. (7), can be naively interpreted in the 
strong nonperturbative way, i.e. X{Q^) ~ l/as{Q^)- Our analysis, however, demonstrates 
that the rise can be explained as ~ InlnQ^, that is natural in perturbative QCD at low 
X (see |]16[, [0-[jl5| and references therein): when the coupling constant is running, the 
renormalization group leads to the small x behavior for the PD ~ ln[as{Q'^)) at the LO 
of perturbation theory and ~ ctsiQ"^) at NLO (see Eqs. (|])-(|1^) and discussions after Eq. 
(0))- 

The good agreement between perturbative QCD and the experiment obtained here and 
in Ref. [0, |l^ demonstrate the fact that for > 2 GeV^ non-perturbative contributions 
as shadowing effects [Q, higher twist effects |^ and others seem to be quite small (see 
also Ref. and references therein) or they seem to be canceled between them and/or 
with ln(l/x) terms contained in the higher orders of perturbation theory. Note, however, 
that higher twist corrections are important at < 1 GeV^, as it has been demonstrated 
in Ref. 



15, 20, 37 



To clear up the correct contributions from nonperturbative dynamics 
and higher orders containing strong ln(l/x) terms, it is necessary further efforts in the 
development of theoretical approaches. 

Moreover, the good agreement between perturbative QCD and experimental data 
at low can be explained with a larger effective scale for the QCD coupling constant 
|T^, [T^. The similar behavior has already been observed in the framework of perturbative 
QCD and in BFKL-motivated approaches [^-[^ (see the recent review in Ref. 



and discussions therein). 



Note that large NLO corrections calculated recently in the framework of BFKL 

lead to a strong suppression of the LO BFKL results for the high energy 

g an H). The careful 



(see also 47 



asymptotic behavior of the cross-section (see, for example, 
inclusion of NLO corrections leads to results which are quite close to those obtained 
in pure perturbative QCD 



This fact can give an additional support for the good 
applicability of perturbation theory in the small x range, where, as was expected before, 
nonperturbative effects should give an essential contribution. 

As the next step, it could be very useful to apply the generalized doubled asymp- 
totic scaling approach to make a combined analysis of HERA data for F2, dF2/dlia{Q'^), 
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din F2/ dln{l/ x) and Fl. We hope to consider this work in a forthcoming article, includ- 
ing higher twist corrections in the Q^-evolution approach given by Eqs. (pD-(pTD. We also 
plan to consider additional terms in the initial condition given by Eq. (|l|) proportional to 
ln(l/x) and ln^(l/x). 

We hope that the analysis will be relevant to find out the kinematical region where the 
well-established perturbative QCD formalism can be safely applied at small x. Moreover, 
the study should clear up the reason of the good agreement between the small x relation 
of Fl, F2 and dF2/dln{Q'^) obtained in pure perturbative QCD in Ref. |^ (and based on 
previous works p6|, ^), the experimental data for these structure functions 0] and 
the predictions of Ref. fsTtl in the framework of fej-factorization 



21, 22, 
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Figure captions 



Figure 1. The derivative (ilnF2/(iln(l/a;) (effective slope A ) as a function of x for 
different values. Data points are from HI Only statistical uncertainties are shown 
The solid line is the result from a fit using in Eq- (15) with fixed Q^ = 1 GeV^ and 



Xo = 1. The dotted line is the same as the solid but with the parameters from a fit to 
F2 data in Ref. |T^. The dashed line corresponds to the asymptotic expression X^p^'"^'^ in 
Eq. dll). 

Figure 2. The derivative (ilnF2/(iln(l/x) (effective slope A ) as a function of Q^. 
Data points are from HI Outer error bars include statistic and systematic added in 
quadrature. Inner bars correspond to statistic errors. The solid line corresponds to the 
HI parameterization given in Eq. (7). Dotted and dashed curves are produced as in 
Fig. 1. For the lower (upper) curves, the value x = 10"'^ {x = 10~^) was used. 

Figure 3. The derivative (ilnF2/(iln(l/a;) (effective slope A ) as a function of Q^. 
Data points are from HI Error bars and solid line are as indicated in Fig. 2. The 
dashed lines were calculated with Eq. ([16|) using x = a ■ 10^^ ■ with a = 0.1, 1 and 10. 
Upper curves correspond to larger x. 
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